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Abstract. In this paper we address the existence, the asymptotic behavior and sta- 
bility in U' and L'''°° , | < p < oo, for solutions to the steady state 3D Navier-Stokes 
I equations with possibly very singular external forces. We show that under certain small- 

ness conditions of the forcing term there exists solutions to the stationary Navier-Stokes 
equations in L'' spaces, and we prove the stability of these solutions. Namely, we prove 
that such small steady state solutions attract time dependent solutions driven by the 
same forcing, no matter how large the initial velocity of non-stationary solutions is. 
We also give non-existence results of stationary solutions in , for 1 < p < |. 

< 

> . 1. Introduction 

In this paper we consider the solutions to the three-dimensional steady state Navier- 
Stokes equations in the whole space M'^, 

(N ' 

00 : (1.1) ' 



(M: • iV-c/ = o. 

^ . Here U = {Ui,U2,Us) is the velocity, P the pressure and / = (/i,/2,/3) a given time 

independent external force. Equation (II. ip will be complemented with a boundary condi- 
■ tion at infinity of the form U{x) — )■ in a weak sense: typically, we express this condition 

0^ . requiring that U belongs to some spaces. Three problems will be addressed. 

^ We will first establish the existence of solutions U £ L^, with | < p < oo, to equations 

(jl.ip for (small) functions / as general as possible, and non-existence results in the range 

X: i<p<i- 

^ ' Next we will study the asymptotic properties as \x\ — )• oo for a relevant subclass of the 

solutions obtained. 

The third problem at hand is the stability of the solutions in the sense of solutions 
to (jl.ip being "fixed point" in LP to the non-stationary incompressible Navier-Stokes 
equations in 

dfU + u ■ Vu + Vp = Au + / 
(1.2) { V-n = 

u{0) = no, 
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where u, p are the time dependent velocity and pressure of the flow. We assume / to 
be constant in time, but our methods could also be applied to the more general case of 
time dependent forces suitably converging to a steady state forcing term. We will show 
that small stationary solutions U of (jl.ip will attract all global non-stationary solutions u 
to ()1.2p verifying mild regularity conditions, and emanating from possibly large data uq. 
This will be achieved by first proving that a wide class of global solutions of ()1.2p must 
become small in L^'°° after some time, and then applying the stability theory of small 
solutions in L^'°° as developed, e.g., in [9l [20l [30]. In addition, for small solutions, we 
will extend the results on the stability in the existing literature by giving necessary and 
sufficient conditions to have — t- [/ in L*' as t — ?• oo. 

The existence and stability of stationary solutions is well understood in the case of 
bounded domains. See for example [10]. For related results in exterior domains we refer 
the reader to [HI [I2l [HI [16] . A wider list of references regarding connected literature can 
be found in [3]. For example, the existence and the stability of stationary solutions in 
with p > n, where n is the dimension of the space, is obtained in [24], under the condition 
that the Reynolds number is sufficiently small, and in |20] . |30j under the assumption that 
the external force is small in a Lorentz space. Similar results in the whole domain M", 
always for p > n, have been obtained also in jl9j . [8], [9]. 

On the other hand, not so much can be found in the literature about the existence 
and stability of stationary solutions in with p < n. This problem have been studied 
recently in the case n = 3 and p = 2 in [3]. In this paper we extend the results of [3] to the 
range | < p < oo, and improve such results also in the case p = 2 by considering a more 
general class of forcing functions. The methods in this paper differ completely from the 
ones used in [3] . In the former paper the construction of solutions with finite energy was 
based on a well known formal observation: if <I> is the fundamental solution for the heat 
equation then ^{t,-) dt is the fundamental solution for Poisson's equation. Using that 
idea it was possible to make a time dependent PDF similar to the Navier-Stokes equation 
with / as initial data with a solution that can be formally integrated in time to find a 
solution of (jl.ip . 

As we shall see, the conditions on / in the present paper which yield that [/ G are, 
essentially, necessary and sufficient. This will be made possible by a systematic use of 
suitable function spaces. 

One could also complement the system (jl.ip with different type of boundary condition 
at infinity. For example, conditions of the form U{x) — )■ Uoo as — )■ oo, where Uoo £ 
and Uoo 7^ are also of interest. However the properties of stationary solutions satisfying 
such condition are already quite well understood. We refer to the treatise of Galdi [H] for 
a comprehensive study of this question. 

On the other hand, the understanding of the problem in the case U^o = is less 
satisfactory. For example, the construction of solutions obeying to the natural energy 
equality (obtained multiplying the equation (jl.ip by U and formally integrating by parts) , 
without putting any smallness assumption on /, is still an open problem. The main 
difficulty, for example when Q = (or when Poincare's inequality is not available), is 
that the usual a priori estimate on the Dirichlet integral 

l|VC/||i2<||/||^_l 
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ensures only that U G C L^: but to give a sense to the integral in the formal equality 

j [V ■{U(g)U)] ■Udx = 

one would need, e.g., that U belongs also to L^. 

More generally, one motivation for developing the theory (especially for low values 
of p) of stationary solutions is that this provides additional information on the asymptotic 
properties of U in the far field. On the other hand, condition like U £ for large p are 
usually easily recovered via the standard regularity theory, as bootstrapping procedures 
show that weak solutions U G are regular if / is so. See also [2H for this case. 

The paper will be organized as follows. After the introduction we have a section of 
general notation, where we recall definitions of several function spaces which will be needed 
in the sequel. 

Section two deals with the existence of solutions in L^', | < p < oo. Section three 
addresses the pointwise behavior in of the solutions and the asymptotic profiles. We 
note that the study of the asymptotic profiles has been largely dealt in the literature, 
starting with the well known results of Finn [13] in exterior domains. Our results being in 
the whole domain are simpler, but we are able to get them with weaker conditions. Non 
existence results of (generic) solutions in [/ G L^, p ^ | will also follow from such analysis. 

Section four handles the stability of stationary solutions. More precisely in the setting 
of the Navier-Stokes equation we investigate the stability of the stationary solution U in 
the and the Lorentz LP'°°-norms. We consider a possibly large non-stationary 
solution and a stationary solution U £ /^3,oo n or [/ G 2^3,oo p ^p,oo -^j^j^;}^ small in 
^3,oo_ show that the non-stationary solution eventually becomes small in (but 
does not converge to in this space), we prove some decay estimates for it and we give a 
necessary and sufficient condition to have that u{t) — )• [/ in or LP'°°. 

The fact that small steady state solutions U attract small non-stationary solutions was 
proved by several authors in different functional settings, see, e.g. [HI O [191 [201 [30] . The 
main novelty of our approach is that we can prove the same result for a class of large 
solutions. At best of our knowledge, this was known only in the particular case U = (see 
[T1[T5]). Our main tool will be a decomposition criterion for functions in Lorentz-spaces. 

1.1. Notations. 

1.1.1. Function spaces. We recall that the fractional Sobolev spaces (or Bessel potential 
spaces) are defined, for s G M and 1 < p < cxo, as 

H; = {fG S\R^) : J--i(l + |^P)t/ G LP}, 

and their homogeneous counterpart is 

H; = {feS'{R^):T-'\^\'feLP}. 

Their differential dimension is s — -. We will only deal with the case s — - < 0, so that 

the elements of Hp can indeed be realized as tempered distributions. As usual, we will 

simply write and instead of i?! ^2 ^'^^ the classical Sobolev spaces. 

The fractional Sobolev spaces can be identified with particular Triebel-Lizorkin spaces, 
namely Fp'^ and Fp'^. This identification will be useful, because it allows us to handle 
the limit case for p = 1: the corresponding spaces are defined as above, but replacing 
with its natural substitute, i.e., the Hardy space Ti^. Similarly, in the limit case p = oo 
one replaces L°° space with BMO. The classical reference for function spaces is |29| . 
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We will make extensive use of the Lorentz spaces L^'^, with 1 < p < oo and 1 < g < oo. 
For completeness we recall their definition. 

Let {X, A) be a measure space. Let / be a scalar-valued A-measurable function and 

Xf{s) = X{x : f{x) > s}. 

Then re-arrangement function /* is defined as usual by: 

rit) = mf{s:Xfis)<t}. 



By definition, for 1 < p < oo, 
LP'''(]R") = {/ : 

where 



C, measurable : 



LP, 9 



< oo}, 



\i q < oo, 



sup{ipr(t)}, ifg 

t>0 



oo. 



We note that it is standard to use the above as a norm even if it does not satisfy the 
triangle inequality since one can find an equivalent norm that makes the space into a 
Banach space. 

In particular, LP'°° agrees with the weak space (or Marcinkiewicz space) 
LP* = {/ : M" ^ C: / measurable, < oo}. 

The quasi-norm 



LP* 



supi[Aj(t)]p 

t>o 



is equivalent to the norm on LP'°°, for 1 < p < oo. 

Our measure A will be chosen to be the Lebesgue measure. The Lebesgue measure of 
a set A will be denoted by mes{A). For basic properties of these spaces useful reference 
are also |31] . |21] . It is well-known that the space L^'"^, 1 < p < oo and 1 < q < oo, is 
the interpolated space L^'*^ = [L^,L°°]^_i . Here [•, „ denotes the interpolated space 

by the real interpolation method. Using the reiteration theorem for interpolation, see |2H 
Theorem 2.2], one has that L^''' = [LPi'^i , L^^'^ajg for all 1 < < p2 < oo, 1 < g, q-2, q-2 < 
oo, < 6* < 1 and i = + A. in particular, one has that n L^^'^a c LP'I for ah 

' p Pl P2 ' 

l<pi<p<P2<oo and 1 < q,q2,q2 < oo. The Holder inequality in Lorentz spaces can 
be stated in the following form. 

Proposition 1.1. Suppose that 

111 111 

1 < p,pi,P2 < oo, I < q,qi,q2 < OO, - = 1 and - = 1 . 

p Pi P2 q qi q2 

Then the pointwise product is a bounded bilinear operator from L^i'"?! x L'p^''^^ to L^''^, from 

lP,q X ^oo ^ J^p,q J f J^p,q ^ J^p' ,q' J^l ^ 1+J. = 1 and ^ + \ = I . 

J p ' p' q q' 

The proof of this proposition can be found in [2H Proposition 2.3]. The similar property 
for convolution is proved in |2H Proposition 2.4] and reads as follows. 



Proposition 1.2. Assume that 

1 < p,pi,P2 < oo, 1 < q,qi,q2 < oo. 



111 ,111 
1 H — = 1 and — = 1 . 

p Pi P2 q qi q2 
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Then the convolution is a bounded bilinear operator from L^i'^i x ^^2,92 L^'i^ from 
LP^i X to LP'i and from 1^^" x L^''^' to L°° where i + A = l and^ + ^ = 1. 

J p ' p' q q' 

We also recall the definition of the Morrey-Campanato spaces. In their homogeneous 
version, for 1 < (7 < p, their elements are all the L]^q^(M'^) functions / satisfying 

\\f\\M, , = sup sup Rv~m dx) ^ < 00 

3;oeM3-R>0 ^J\x-xo\<R ' 

We recall that 

(1.3) LP = LP'P = Mp^p C C A^p,5, 1 < g < p < 00, 

with continuous injections. The M.p^q spaces are of course increasing in the sense of the 
inclusion as q decreases. On the other hand, the L^'"?- spaces increase with q. 

For > 0, we introduce the space Eq of all measurable functions (or vector field) / in 
M^, such that 

ll/lUa = esssup^gjj3 \xf\f{x)\ < 00. 

1.1.2. Other notations. We denote by P = Id — VA^^div the Leray projector onto the 
divergence-free vector field. Notice that P is a pseudodifferential operator of order zero, 
which is bounded in Hp, Hp and L^''^, for l<p<oo,l<g<oo and s G M. Thus, 
when / belongs to those spaces, the validity of an Helmholtz decomposition / = P/ + Vg 
implies that one could assume, without restriction, that / is divergence-free. 

However, we will not make this assumption in order to avoid unpleasant restrictions, 
especially when working in weighted spaces (notice that P is not bounded in Eq) or in L^. 
Indeed, it has some interest to consider integrable external forces with non-zero mean, 
which prevents div / = 0. 

2. Solutions in Lp{R^) 

The equations (jl.ip are invariant by the natural scaling {U,p,f) 1— )• {U\,px, f\) for all 
A > and C/a = AC/(A-), -Pa = A2p(A-) and fx = \^f{\-)- Following a well established 
procedure, not only for Navier-Stokes, we consider the following program: 

(1) Existence: first construct (rough) solutions [/ in a scaling invariant setting, i.e. in 
a functional space with the same homogeneity of assuming that the norm of / is small 
in a function space (as large as possible) with the same homogeneity of . 

(2) Propagation: deduce from additional properties of / (oscillations, localization,. . .) 
additional properties for U (localization, asymptotic properties,. . . ). 

We will not discuss the propagation of the regularity since this issue is already well 
understood (see |14]). For example for, not necessarily small, external forces belonging to 
H~^ n -ff*, with s > |, one deduces that solutions with finite Dirichlet integral are twice 
continuously differentiable and solve (jl.ip in the classical sense. 

Concerning the first part of this program, in order to give a sense to the nonlinearity 
one wants to have U E Lfoc- noticed in [22], the largest Banach space X, continuously 
included in L^q^(]R^), which is invariant under translations and such that ||f/A||x = ||t^||X)is 
the Morrey-Campanato space A^3,2- Therefore, the weakest possible smallness assumption 
under which one can hope to apply the first part of the program should be 

l|A-VlU3,2 < £■ 
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However, it seems impossible to prove the existence of a solution under this type of 

condition. Indeed, U ®U would belong to i, and the singular integrals involved in 

2 ' 

equivalent formulations of (jl.ip are badly behaved in Morrey spaces of L\^^ functions (see 
the analysis of Taylor [27j and in particular Eq. (3.37) of his paper). 

Here the situation is less favorable than for the free non- stationary Navier-Stokes equa- 
tions, where the existence of a global in time solution can be ensured if the initial datum of 
the Cauchy problem is small in A^3,2 (or even under more general smallness assumptions, 
see [2I])- The complication, in our case, arises from the lack of the regularizing effect of 
the heat kernel. 

On the other hand, the above difficulty disappears in the slightly smaller spaces A^3,g- 
Indeed, Kozono and Yamazaki established the following result 

Theorem 2.1 (See [19]). Let 2 < q < 2>. Then there exists a positive number 5q and a 
strictly monotone Junction 0Jq{5) on [0,(5^] satisfying ojq{Q) = 0, such that the following 
holds: 

• For every f G 'D'(M^) there exists at most one solution U in Mz,q satisfying 

• For every tempered distribution f such that A~^/ € Als.g, o-nd 5 = || A^"'^/||_a43 ? 
6q, there exists a solution U G A^3,g of (jl.ip . such that \\U\\M3q < ^q{S)- 

This result provides a satisfactory answer to Part 1 of the above program, but it seems 
difficult to make progress in Part 2 using such functional setting. For example, a very 
strong additional condition like / € 5o(M^) (the space of functions in the Schwartz class 
with vanishing moments of all order), and / small, but only in the A^s^g-norm (with 2 < 
q < 3), seems to imply no interesting asymptotic properties for U (such as U £ L^ with 
low p). 

On the other hand the ^As^q spaces, as g t 3, become very close to L^'°° as can be 
seen from relation ()1.3p . The purpose of our first theorem is to show that one can obtain 
propagation results according to Part 2 of our program, by strengthening a little the 
smallness assumption, and requiring that 

(2.1) ||A-V||i3,^ <ei. 

The continuous embedding of L^ into the weak space L^'°° implies that condition (j2.ip 
will be fulfilled if, i.e., f G H^^ with small H^^-norm. Moreover, the continuous embed- 
ding 

3 • o 

3 . 3 

shows that the case of forces / E H^2 with small i7^2-norm is also encompassed by (|2.ip . 
We now state our first theorem. 

Theorem 2.2. There exists an absolute constant ei > with the following properties: 

• If f £ 5'(M^) is such that A~^f € L^'°° and satisfying condition (12. ip . then there 
exists a solution U € L'^'°° of (jl.ip such that 

(2.2) ||f^||L3.- < 2||A-ip/||i3,oo. 

(The uniqueness holds in the more general setting of Theorem \2.1\) . 

• Let I < p < 00. If U is the above solution then we have more precisely 

U G n LP if and only if Ff G Hp'^. 
In this case (and if p ^ 3) U £ L'^ for all q such that 3 < q < p (or p < q < 3). 
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Moreover, U belongs to L'^'°^ n (respectively, U £ L^'°° n BMOJ if and only 
ifA-^Ff G (respectively, A'^Ff G BMO). 

Remark 2.3. Important examples of solutions that can be obtained through this theorem 
are those corresponding to external forces / = (/i, /2, /s) with components of the form e6, 
where 5 is the Dirac mass at the origin. Notice that / H~^/^. However, assumption (j2.ip 
is fulfilled, because A~-'^/(x) = jfj (ci, C2, C3). 

In fact, due to the invariance under rotations of (jl.ip . in this case one can always fix 
a coordinate system in a way such that / = (eJ, 0,0). The solutions that one obtains 
in this way are well-known: they are the axi-symmetric solutions (around the xi axis) 
discovered by Landau sixty years ago, with ordinary differential equations methods. These 
are solutions that are singular at the origin — in fact the components of the velocity field 
are homogeneous functions of degree —1 — and smooth outside zero. They can also be 
seen as self-similar stationary solutions of the non- stationary Navier-Stokes equations. 

We refer to [S] for an explicit expressions and other interesting properties about these 
solutions and to [26] (see also [28]) for related uniqueness results. 

Remark 2.4. The particular case p = 2 is physically relevant since it corresponds to finite 
energy solutions. The conclusion U £ was obtained by Bjorland and Schonbek [3], un- 
der a technical smallness assumption non invariant under scaling. Part (2) of Theorem l2.2l 
improves their result. Indeed the same conclusion can be reached under the more general 
conditions (12. ip and / G H~^. In particular, it follows that / G i?~2 n with / small 
in H 2 would be enough to get U G L . This fact was pointed out to the first and the 
last author by an anonymous referee of their paper [3]. 

Roughly speaking, for / G H~2 ^ the additional requirement / G (which turns out 
to be also necessary for obtaining U G L^, up to a modification of / with an additive 
potential force, which in any case would change only the pressure of the flow), is formally 
equivalent to the additional vanishing condition /(^) = o(|^|2) as |.^| — )■ 0. 

Remark 2.5. The first conclusion of Theorem 12.21 bears some relations with the work 
of Kozono and Yamazaki [20] and Yamazaki [30], where they also obtained existence 
results of (possibly non-stationary) solutions in Lorentz-spaces and in unbounded domains. 
However, the assumptions in [201 [30] on the external force reads / = div F, with F small 

3 

in L2'°°. Their condition is more stringent than our condition (j2.ip because it involves 
more regularity (one more derivative, or more precisely, one less anti-derivative) on /. 

The first part of Theorem 12.21 is also related to the work of Cannone and Karch [9]. 
There, the authors constructed non-stationary solutions of Navier-Stokes in the whole 
space in L^(L^'°°) with initial data small in and external force such that 

f e^^-'^^¥f{s)ds 

Jo L3.°o 

is small. With some modifications of their proofs it would be possible to deduce the first 
conclusion of our theorem from their result, by considering time-independent external 
forces (in this case the above condition boils down to (j2.ip ). We prefer however to give 
a self-contained proof directly in the stationary case, because this allows us to obtain 
necessary and sufficient conditions. Moreover, none of the these papers addressed the 
construction of solution in with p < 3. 

We recall a well known fixed point Lemma for bilinear forms that will be needed in the 
sequel. The proof can be found in 
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Lemma 2.6. Let X be a Banach space and B : X x X ^ X a bilinear map. Let \\ ■ \\x 

denote the norm in X. If for all xi,X2 £ X one has 

\\B{xi,X2)\\x < 'n\\xi\\x\\x2\\x- 

Then for for all y £ X satisfying 47/||y||x < 1, the equation 

X = y + B{x, x), 

has a solution x € X satisfying and uniquely defined by the condition 

\\x\\x < 2||y||x- 

Remark 2.7. The proof of this lemma also shows that x = lim where the approximate 

solutions Xk are defined hy xq = y and = y + i?(xfc„i, Xfc„i). Moreover < 2||y||x 

for all k. 

Proof of Theorem \2.2[ 

We use a method of mixed bilinear estimates, inspired from |17j . Let us set 
Uo = -A-^Ff, B{U,V) = A-'^FV ■ {U (g)V). 
Then the system (jl.ip can be rewritten as 

(2.3) U = Uo + B{U, U) 

and the solutions of this equations are indeed weak solutions of (jl.ip . This equation can 
be solved applying the standard fixed point method as described in Lemma 12.61 in space 
We have the estimate 

(2.4) \\B{U,V)\\iji,oo < (:7i||C/||i3,oo||T/||i3,oc, 

for some Ci > independent on U and V. Note that an estimate similar to ()2.4p has been 
proved e.g. by Meyer in [22j in the case of the non-stationary Navier-Stokes equations 
(the bilinear operator B is slightly different in that case). 

To prove (j2.4p . we only have to observe that the symbol m-(^) of the pseudo-differential 
operator A~^Pdiv is a homogeneous function of degree —1, such that m(^) G C°°(M^\{0}). 
Thus, the corresponding kernel m is a homogeneous function of degree —2, smooth outside 
the origin (more precisely m — (j^j,h,k}j,h,k=i,2,3 

and mj h^k are homogeneous functions of 

degree —2). 
In particular, 

B{U,V) = m{D){U ®V) with m G Li'"" 

Thence, 

(2.5) \\m{D)v\\LV2...<C{pi,qi)\\v\\Lvi^.,, ^ = J- + | _ 1, 

by the Young inequality stated in Proposition 11.21 Applying this to v = U ®V and 
using Proposition 11.11 to deduce that for U,V ^ ^3,oo ^i^s that v G L2'°°, we get 
estimate (12. 4p with Ci = C(|,oo). Hence by Lemma 12.61 it follows that, provided that 
4||C/o||Ci < 1, there exists a solution of ()2.3p satisfying ()2.2p . 

To prove Part 2, we make use of approximate solutions of ^(U) = Uq + B{U, U). That 
is we choose a sequence satisfying Uk = Uq + B{Uk-i, Uk~i) and use Remark 12.71 to state 
that C/fe — [/ in L^'°° as /c — )• oo. We show now that 

(2.6) \\B{UkM\\Lv<C{p)\\Uk\\v^,^\\Uk\\Lv, |<p<oo. 



K Pi < 3, 
1 < (7i < oo. 
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valid for some positive function p i— )• C{p), continuous on (|,oo). To obtain this estimate 
we use the Holder inequality given in Proposition 11.11 to deduce that (8) C/fc|| sp < 

C2(p)||?7fc||Lp||f^fc||L3.°°- Relation (|2.5p for pi = ^2- and qi = p completes the proof of 
(ESI). 

By Part 1, applied to the approximations Uk, we know that ||{7fc||^3,oo < 2||J7o||l3,oo. 
Choose ||C/o||l3,cxj < coei, then we get from (j2.6p . for | < p < oo 

(2.7) \\Uk+i\\Lp < \\Uo\\Lp + 2coC{p)ei\\Uk\\Lp. 

If Ff G Hp^, then Uq G and so, by induction, ||i7fc||LP < oo for all k. Provided 
2coC{p)£i < 1, iterating inequality (j2.7p implies that Uk is uniformly bounded in with 
respect to k, and hence U S L^. 

However, C{p) blows up as p — ?• | or p — )• oo, and we want to have a smallness assump- 
tion independent of p. To circumvent this difficulty, we replace, if necessary, the constant 
£i of Part 1 of the theorem with a smaller absolute constant (still denoted £i), in a such 
way that 2coei < 1/ sup2<p<7 C{p). Then the above argument yields the conclusion of the 
"if part" of the theorem in the case 2 < p < 7. To prove the "only if" part one simply 
uses estimate (|2.6p with Uk = U together with 12. 3p to get Uq = — A(P/) e L'P, hence 

Let us now consider the case Ff G H"^, | < p < 2. Then [/q G n L^'°° and 
by interpolation Uq G L , so using the case 2 < p < 7 we get that U £ L^. On the 

3 

other hand, according to Proposition 11.21 the space L2'°° is stable under convolution with 
L^-functions so 

B{U, U) = m{D){U [/) G Lt'°°. 

But from estimate ( 12.41) we know that B{U,U) G L^^°°. By interpolation, B{U,U) G LP. 
Combining this with equality (j2.3p yields U G L^. Conversely, suppose that U G L^. Since 
we already know that the solution U G -L'^'°°, we deduce by interpolation that U £ L'^. 
The argument above shows that B{U,U) G L^. Hence by (|2.4p it follows that Uq G L^, 
and this, in turn is equivalent to Ff G H'^. 

We now consider the case Uq G with 7 < p < oo (this is equivalent to Ff G 
if p < oo). Since Uq G L^'°°, we have by interpolation that Uq L'^ (1 L''. From the 
previous case, we infer that U G n L'^. By interpolation, we also have that U G L^'^. 
Using Proposition 11.11 this implies that U^U G L^'^, so from Proposition II . 21 and recalling 

3 

m G L2'°° we get that 

B{U, U) = m{D){U ®U)eL°°. 

But we also know that B{U, U) G so by interpolation B{U, U) G . From (fO]) we 

conclude that U G L^. The same argument also shows that U £ implies Uq G L^. 

Finally, the BMO case follows in the same way. Indeed, the argument above shows that 
if Uo or U belong to BMO, then B{U,U) G But C BMO, so B{U,U) G BMO. 
From relation (j2.3p we see that U G BMO iff Uq G BMO. This completes the proof of 
Theorem [221 □ 



Remark 2.8. With the same proof, one can show the following equivalent condition for 
the stationary solution U constructed in Theorem 12.21 to belong to L^'^. If p G (|, oo) and 
r G [1, oo] then U G LP''' if and only if A'^Ff G LP'' 
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3. Pointwise behavior in and asymptotic profiles 

In the previous section we dealt with forces such that A~^/ G Since the typical 

example of a function in L^'°° is it is natural to ask which supplementary properties 

are satisfied by the solution when |A^^/(x)| < The theorem below provides a 

rather complete answer. 

In particular, we will obtain exact asymptotic profiles in the far field for decaying 
solutions of (jl.ip . Starting with the work of Finn (see [13] and the references therein), a 
lot is known about the spatial asymptotics of stationary solutions in unbounded domains. 
The case of the whole space that we treat in this section is of course simpler than the case 
of exterior domains or aperture domains considered e.g. in (14j . Nevertheless, focusing on 
this case allow us to put weaker (and more natural) smallness assumptions on the force, 
thus providing a more transparent presentation of the problem. 

We observe here that, despite the unboundedness of P in the Eq spaces, it is fairly easy 
to ensure e.g. that A"ip/ G El. Indeed, one has for example that 

(3.1) ||A-ip/||^^<C(||/||^3 + 11/11^0- 

Notice that all the norms in inequality (j3.ip are invariant under scaling. The above 
inequality can be proved with a simple size estimate (using that A^^P is a convolution 
operator with a kernel m satisfying |wi(x)| < C|x|~^). The same conclusion A~-'^P/ G Ei 
can be obtained also via the Fourier transform (using classical results in [25]), assuming, 
e.g., f = V ■ F where -F is a two dimensional tensor with homogeneous components of 
degree —2, smooth outside the origin. 

Let us recall the imbedding Ei ^ L'^'°°, thus the smallness assumption in the space Ei 
implies a smallness assumption in L^'°° . 

The spirit of Theorem 3.1 below is close to a previous work of the second author (see [6]) 
in which similar conclusions are shown for the time-dependant Navier-Stokes equation in 
the whole space. 

Theorem 3.1. There exists an absolute constant £2 > (with £2 a priori smaller than 
the constant ei of Theorem \2.2\) such that: 

• If f € S'iR-^) is such that A'^P/ G Ei and ||A-ip/||^^ < 82, then the solution 
U G L^'^ obtained in Theorem \2.S\ satisfies 

ni^^<2||A-ip/ii^^. 

• Let < 9 < 2. Under the additional assumption A~"'^P/ G Eg, we have also 
U €Eg. 

• In particular, if A~^P/ £ Eq Ci E2, with small Ei-norm, then U satisfies the 
pointwise estimate 

\U{x)\ < C(l + |x|)-2. 
In this case the solution U has the following profile as \x\ — t- 00; 

(3.2) U{x) = -A-'^¥f{x) + m{x) : (^j U®U^ + 0{\x\-^ log(|x|)) , 

where m = {mj^h,k) is the kernel of A~^F div and mj^h,k{x) c^'^e homogeneous func- 
tions of degree —2, outside zero. Furthermore, 

(3.3) m{x) : (^J U =0 if and only if 3 c G M s.t. J UhUk = c5h,k, 
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for h,k = 1, 2, 3, where 5^^^ = or 1 if h ^ k or h = k. 
Remark 3.2. Let us be more explicit with our notation: by definition, for j = 1, 2, 3, 

7n(x): j\u®U)^^ = ^ mj^hAx)(^j Uh{y)Uk{y)dy^. 

Moreover rnj^fi,ki^) = dh^ii,k{^)^ where Mj^^ is the tensor appearing in the fundamental 
solution of the Stokes equation. The computation of M goes back to Lorentz (1896). 
See \\.4\ Vol. I, p. 190] for the explicit formula. 

Remark 3.3. For example, it follows from this theorem that, if / € is such that 

supp/ and / satisfies the previous smallness assumption, then 

U [x] ~ m{x) : U <^U^ , as|x|— s-oo 

provided that the right-hand side does not vanish. Indeed, we have in this case A~^P/ G 
5(M^). In particular |f/(x)| < C(l + But the improved estimate U{x) = o(|x|"^) 

as \x\ —7- oo holds if and only if the flow satisfies the orthogonality relations (13. 3p . Of 
course, generically it is not the case. This implies the optimality of the restriction < 2 in 
Theorem 13.11 as well as the optimality of the restriction p > | appearing in Theorem 12. 2[ 

It is possible to relax the condition that supp/ assuming, instead that |/(^)| < C\£,\^ 
for a sufficiently large A: > 0. As noticed in [3], this is essentially an oscillatory condition 
on /, describing the large time behavior of the solution of the Cauchy problem for the 
heat equation. 

Remark 3.4. Examples of (exceptional) stationary flows satisfying the orthogonality re- 
lations (|3.3p . and such that U{x) = 0(|a;|~^ log(|x|)) , are easily constructed by taking 
/ satisfying the assumptions of the previous remark and additional suitable symmetries. 
An axi-symmetry condition would not be enough: one rather needs here polyhedral- type 
symmetries. The suitable symmetries to be imposed on / can be classified exactly as done 
in [Ij, in the case of the non-stationary Navier-Stokes equations. For example the two 
conditions Rf{x) = f{Rx) and Sf{x) = f{Sx) where R, S are the orthogonal transforma- 
tions in M.^ R: {xi,X2,X3) ^ (x2,X3,xi) and S: {xi,X2,X3) i— {—xi,X2,X3) are sufficient. 
See [1| for explicit examples of this type of vector fields. 

On the other hand, explicit examples of solutions U = Uj which do not satisfy the 
orthogonality relations can be obtained simply by taking / = ry/o with rj > sufficiently 
small and /o S satisfying the conditions of Part 2 of Theorem 12.21 with p = 2 (this 
implies that [//q G L^). If, in addition, there is no c (zM. such that 



J {A~^Ffo)h{A-^rfo)k = c6h,k 



then Uf cannot satisfy the orthogonality relations, provided > is small enough. The 
proof of this claim relies on an argument that has been used in [5] in the setting of the non- 
stationary Navier-Stokes equations. These observations lead us to the following theorem, 
containing the announced non-existence result of generic solutions in L^, p < |. 
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Theorem 3.5. Let fo = (/i,/2,/3) be a divergence- free vector field such that f G 
and supp{f). Assume also that the matrix 

ifo)jifo)k ,A 

is not a scalar multiple of the identity. Then there exists r]Q > such that the solution 
of p.ip with f = 77/0 and < 77 < 770 satisfies, 

(3.4) c\x\~^ <\U{x)\ <C\x\-^, |x| » 1, 

where C > is independent on x and c = c(|^) is independent on \x\; moreover, c(|f|) > 

on a set of positive surface measure on the unit sphere. In particular, U LP(]R^) for all 
3 
2 ■ 



1 < p < 3 



Proof of Theorem \3. 1[ We already have, by Theorem 12.21 a solution in To see that 

such solution belongs more precisely, to Ei we only have to prove the estimate 

(3.5) \\BiU,V)\\^^<C\\U\\^jV\\^^, 

for some C > independent on U and V. Indeed, an application of Lemma 12.61 shows 
the existence and the uniqueness of the solution C7 in £^1. This solution also belongs to 
since Ei C L^'°°. Of course, the re-application of the fixed point argument requires 
that we replace the constant ei > of Theorem 12.21 by a smaller one. Relation (j3.5p is a 
particular case of the following lemma: 

Lemma 3.6. Let 61,62 be two real numbers such that 1 < 9i + 62 < 3. There exists a 
constant C such that 

Moreover 

\\B{U,U)\\^^<C{\\U\\%^+\\U\\l). 

Proof. Recall that B{U,V) = m * (U V) with m homogeneous of degree —2. Since 
|?72(x)| < C|x|~^ we can bound 

V) = j mix -y):{U^ V){y) dy < C\\U\\^^^ \\V\\^^^ j _ y\\y\e,^e.. dy. 

It is easy to show that the last integral is a function of | x | homogeneous of order I — 61 — 62, 
so it can be bounded by C\x\^~^^~^^ . 
To show the second part, we decompose 

B{U,V) = ([ +[ )m{x-y):{U®V){y)dy = h+l2. 

We have 

where we used the same scaling argument as above to deduce the last inequality. Next, 
we write for Ii 

\h\<C f —^\U{y)\Uy<^f \U{y)\''dy<^\\U\\l,. 

□ 
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Let us now prove Part 2 of Theorem 13. li We have the additional information A^^Ff £ 
Eq. We argue as in the proof of Theorem 12. 2[ That is we define = [/q + B{U, U) and 

we choose a sequence satisfying = ^{Uk-i)- From Lemma 13.61 we have the estimate 

\\B{Uk, Uk)\\E, < Cgm\\E,\m\^^, 0<e<2, 

for some positive function 9 i— )• Ce, continuous on (0,2). As in Theorem 12.21 part 2 it 
fohows that the sequence of approximate solutions Uk remains bounded in Eq, provided 
that A-ip/ e Ee, for some 6 G (0, 2), and 

2C,||A-ip/||^^<L 

The continuity of Cq allows to obtain the conclusion of the theorem (with a smallness as- 
sumption independent on 9), at least for e.g. 9 G [^) |]- We had to exclude a neighborhood 
of 9 = and of 9 = 2, where Ce blows- up. 

In the case | < < 2, we know that A~^P/ G EiDEq C EiDE?. So, from the previous 

case we deduce that the solution U satisfies U £ EiCi Et_ C L'^ Ci E3 . Using again Lemma 

4 2 

ESI we infer that B{U, U) G E2- But we also know that B{U, U) G Ei so B{U, U) G Eg. 
The conclusion in the case j < 9 < 2 now follows from equation (j2.3p . 

It remains to consider the case 0<9 <l (the case 9 = is contained in Theorem 12.21 

since Eq = L°°). As above, we show that U G Ei D Ei so U £ Ee+i . From Lemma 13.61 we 

2 2 

get that B{U, U) £ Eg so U £ Eg. The proof of Part 2 of Theorem 13.11 is now completed. 

Let us prove Part 3. We will show using decay properties of m and a Taylor expansion 
that for any solution such that |C/(x)| < C(l + we have 

(3.6) A-^FV ■ (U <E,U){x) = m{x): j U <E, U + O {\x\-^ log{\x\)) , as |x| ^ 00. 
But, 

A^^PV • {U (g) U){x) = [ m{x -y):U® U{y) dy 



=m(x): I U ®U - m{x): I U ®U 

J\y\>\A/^ 



+ / [m{x — y) — ra{x)] : U ®U{y)dy 

'l<kl/2 



+ / m{x — y): U ®U{y) dy 

'\x-y\<\x\/2 



+ / m{x — y): U ®U{y) dy. 

J\y\>\x\/2, \x-y\>\x\/2 

The only properties on the kernel m that we will use are |m(x)| < C|x|~^ and |Vm(x)| < 
C|a;|~^. We need to show that all the terms on the RHS of the last inequality (excepted the 
first one) are bounded by C|2;|~^ log for large This follow easily since U £ LP' r\ E2. 
For large the second, the fourth and the last term on the right-hand side are in fact 
bounded by C|x|~^. The third term is bounded by C|3;|~^ log for large \x\, as it can be 
checked applying the Taylor formula to m. This implies both the asymptotic profiles (j3.6p 
and (1321) 
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To conclude, it remains to show that the homogeneous functions 

^ruj^hM^) UhUk, j = 1,2,3, 
h,k 

vanish identically if and only if the matrix J U U is a scalar multiple of the identity. 
We reproduce a computation similar to that in [6l|23]: taking the Fourier transform, the 
above vanishing condition is proved to be equivalent to 

^JAkiO I UhUk = Y.]^ I ^'J^h - Yl / UhUk = 0, for a.e. ^ G K^. 



h,k h h,k 

The conclusion is now obvious. 



□ 



We end this section with the proof of Theorem 13. 5[ 
Proof of Theorem 13.51 We start by choosing 770 sufficiently small such that 

r/o||A-ip/o||i3,^ < ei, 

where ei is the smallness constant from Theorem 12.21 According to Theorem 12. 2|, for 
< ?7 < 770 there exists a unique solution U G 2^3, 00 p ^2 q£ (jup with / = r//o such that 
\\U\\ ^3,00 < 2r/||A ^P/o|| j;^3,oo . It suffices to show that the orthogonality relations (|3.3p does 
not hold true for U. 

Let Wo = — A~^P/o and Uq = rjWo- The hypothesis implies that the matrix JWo^Wq 
is not a scalar multiple of the identity. This means that there exists j ^ k such that either 
JW^Wq 7^ or J|Wjp / /l^o'Pi where denotes the j-th component of Wq. We 
will suppose that / ^0^0 7^ other case being entirely similar. 

We have 

(3.7) I J UW^ - j Uiu^\ = I j {W - ui)U^ + j ui{U^ - U^) 

< ||[/-[/o||i2(||[/||i2 + ||C/o||i2) 

From ()2.3p and (|2.6p with p = 2 and Uk replaced by U we deduce that 

(3.8) \\U-U4l2 = \\B{U,U)\\l2<C{2)\\U\\l2\\U\\l^,o. < 2C(2)r?||H^o||L3.- 
Therefore 

\\U\\l^ < ||C/o||L2 + ||C/-f/o|U2 <7?||M^o||L2+2C(2)r/o||Tyo|U3,.o||f/||i2. 

If we further strengthen the smallness assumption on rjQ by 

1 



< 



4C(2)||H/o||l3,- 
we get that 

\\Uh2<2rj\\Wo\\L^. 
Relation (j3.8|) combined with the previous estimate implies that 

\\U - Uoh^ < AC{2)ij^\\Wo\\l^,^\\Wo\\l^. 
Using the two previous bounds in (j3.7p implies that 



<12C(2)r/=^||I^o||L3-||VFo||i2. 
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Finally 



> V 



l2C{2)rf\\Wo\\L^,o.\\Wo\\l2>^ 



if we further assume that 



m < 



l/^J^ol 



24C(2)||T^o||L3,.o||T^o|li2 



□ 



4. Stability of the Stationary Solutions 

Consider now a mild formulation of the Navier-Stokes equations with time independent 
forcing function / satisfying, as usual, to a smallness condition as in (|2.ip . 

(4.1) «(t) = e*^no+ f e^^~'^^¥fds- f e'^^-'^^FV ■ {u ® u){s) ds. 

Jo Jo 

The two main goals of this section are the following. First we want to establish conditions 
on uq to ensure that the above system has a solution u £ L°°(]R_|_, L^'°°). Next, we want 
to find the largest possible class of solutions u to (j4.ip for which we can say that u{t) 
converges to the steady solution U given by Theorem 12.21 corresponding to the same 
force /. This class will be general enough to include non-stationary solutions in L^'°° 
with large initial data. We will show in particular that a priori global solutions, verifying 
a mild regularity condition but initially large in L^'°°, become small in after some 

time. Only the singularity at infinity of the initial velocity needs to be small in some sense 
which is made rigorous in (j4.17p . For example, we allow an initial velocity uq bounded by 
C/|a;| everywhere and bounded by e/|x| for large x, with C arbitrary and e small. 

We recall that a priori large non-stationary solutions in Bp^q ^ and VMO of the 
Navier-Stokes equations without forcing are known to converge to zero in these spaces 
(see |15l [T]). However, in our case, convergence to zero will not necessarily hold true for 
\\u{t)-U\\ ^s.tx), due to the fact that the smooth function in 5(M^) are not dense in L^'°°. 
Thus, only weaker convergence results should be expected. 

Theorem 14 . 3 1 collects our results on the stability of small solutions u, extending, for flows 
in with time independent forcing term, those of O El [2U1 130] to the case | < p < 3, 
and providing some additional information also for p > 3. Theorem 14.71 contains the 
convergence result of large solutions u to small stationary solutions U. Its proof relies 
on some energy estimates inspired by [15^ [1] and on the results on the stability of small 
solutions prepared in Theorem 14.31 

To begin we first recall a lemma which will be useful for estimating the integral terms 
on the RHS of (]4.ip in spaces. We notice that for the case p = 3, q = 3/2 the Lemma 
below was obtained in several papers, among them the first seems to be in Yamazaki's 
paper |30j . Variants of this lemma can also be found in |22j . in a slightly less general form, 
and in |21j . 

Lemma 4.1. Given any p £ (|,oo) let q = For < a < t, the operator 

L^{(j)){t) = f e(*-^)^PV • (j){s) ds 
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satisfies 

(4.2) \\L^mt)\\LP,^ < C{p) sup ||<^(s)||l,,.o 

0<s<t 

where C{p) denotes a constant independent of a. 

Proof. Let F{t) be the kernel of the operator e*'^Pdiv. First recall the rescaling relation 

F{x,t) = t~^F{x/Vi,l) 

and that E L^nL"^. 

We consider separately the following two pieces. 

f-t rt-X* 

Ai= F{t - s)* (j){s) ds and A2 = F{t - s) * (j){s) ds. 

Jt-X* Ja 

The idea of the estimate is to find, given any fixed A, a A* so that |yl2| < A/2. With this 
choice of A* we can estimate the Lebesgue measure of the set {x : io-(</')| > A} in terms 
of Ai only. In that direction we establish two preliminary estimates. The first is a an 
application of Young's inequality stated in Proposition 11.21 

i-t-X' 

P2IIL- < C / \\F{t-s)\\Lc.4cl)\\L^,oods 



Here, a = 5^5^. The estimate \\F{t — s)\\];^a,i < C{t — s) (that follows from the 

rescaling properties of F) implies 

(4.3) <C(p)(A*)-i^||0||^. ■ 



^9 



Here we have introduced the notation Xq' = L°°{(a,t), L'^'°°). Similarly, \\F(t — s)\\ii < 
{t — s)~2 and 

(4.4) \\Ai\\l.,^< [' \\F{t-s)\\LimL^,^ds<{X*)'^\ 



Jt-X* 

We proceed with the bound for ||Lo-((/>)||lp'°°. Using the definition of the norm and the 
triangle inequality, 

1 

||Lcr(0)(t)||iP,oo < supAmes{x : \Ai\ + \A2\ > A}p 

A>0 

For each A > we may choose A* such that the RHS of (|4.3p is equal to A/2. With this 
choice of A*, 

1 

Ames{a; : \Ai\ + \A2\ > A}p < Ames{a:; : |^i| > A/2}p 

Also, using (jO|) : 

Ames{x : \Ai\ > X/2}p < X'^'i\\Ai\\l,,^ < C||</.||^.,t^ 



Taking the supremum over all A > establishes ()4.2p . 

□ 

The following lemma concerns the large time behavior in L^'°^ of solutions of the heat 
equation. It will provide a better understanding of the statements of our two next theo- 
rems. 



Lemma 4.2. Let f ^ L 



3,00 



LP-SOLUTIONS OF THE STEADY-STATE NAVIER-STOKES EQUATIONS (June 1, 2010) 17 

• Let e > be arbitrary. Then f can be decomposed as / = /i + /2 with fi G and 
||/2||l3,oo < e if and only i/limsupi?mes{|/| > R}^ < e. 

• // Ui-n i?mes{|/| > R}^ = then e^^f ^ in L^^°° ast^oo. 

• There exists some g € L^'°° such that e^^g — )■ in L'^'°° as i —>• oo and such that 
limsupi?mes{|(7| > i?}3 / 0. 

Proof. Assume first that / = /i + /2 with fi G and ||/2||l3,oo < e. We estimate 

uiesm>R}<^ [ 
R Jr3 

so that lim sup i? mes{ I /i I > i?}^ = 0. We also have that 

R-^o 

hmsupi?mes{|/2| > R}^ < sup mes{ | /2 1 > R}^ = ||/2||l3.oo < e. 
R-^o R>o 

Let S G (0, 1). Since {|/| > R} C {|/i| > 6R} U {I/2I > (1 - 6)R} we have that 
limsupi?mes{|/| > R}^ < limsupi?(mes{|/i| xJi?} + mes{ | /2 1 > {1 - 6)R})^ 



R^O R-^O 



< Umsupi?mes{|/i| > (^ii}3 + hm sup ii mes{ | /2 1 > (1 - 6)R}3 
i?-s>0 R^O 

= - lim sup i? mes{ I /i I > ii}3 -| lim sup mes{ | /2 1 > R}^- 

R^O 1 — R^o 

Leting (5—7-0 imphes that limsupi?mes{|/| > R}s < ||/2||i3,oo < e. 

R-^o 

Conversely, assume that limsupi?mes{|/| > R}^ < e. There exists R^ such that 

R^o 

sup i?mes{|/| > R}^ < e. 

0<R<Re 

We set /i = fX{\f\>R^} and /2 = fX{\f\<Re} where x denotes the characteristic function. 
Clearly \f2\ < Re a^d\f2\ < \ f\ so that ~ 

11/211^3,00 = sup i2mes{|/2| > R}^ < sup i?mes{|/| > R}^ < e. 

0<R<Re 0<R<Re 



It remains to show that /i G L2(m3). Let eZhe such that R^ > 2^^. Then 

{\f\>Re}C U {2" < I/I < 2"+!} 



n>N. 



so 



p p 00 « 00 

/ \fi\'= l/l'<E/ |/|2< 4"+imes{2-<|/|} 



< E ^ii/iii^-<°°- 

This shows the first part of the lemma. 



00 . 
4 



2*^ 

n=N^ 
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Assume now that Mm i?mes{|/| > R}^ = and let e > be arbitrary. Using the first 

part we decompose / = /i + /2 with fi E and ||/2||l3,oo < e. The standard decay 
estimates for the heat equation implies that ||e*'^/i| 1^,3,00 < Ct~3 ||/;^||^2 — )• as t — )■ oo. 
Moreover, ||e*^/2||L3,oo < ||/2||l3,oo < e. We infer that 

limsup ||e*^/||i3,oo < limsup(||e*^/i||/^3,(x) + ||e*^/2||i3,(x>) < e. 
Letting e — )• yields limsup ||e*'^/||j;^3,oo = 0, as required. 

t—>oo 

To prove the third part of the lemma, we choose 

9ix) = ^^, (x) = {1 + \x\')2 . 

[X) 

It is a straightforward calculation to check that 

1 / 47r . 1 

limsup iimes{|5(| > Rjs = (— j^. 

On the other hand, we will show that e^^g G which by the decay estimates for the 
heat equation implies that ||e*^(7||i3,oo < C{t — \)~*\\e^^g\\L^ — )• as t — )• oo. 
Since the kernel of the operator e*^ is 7r~2e~l^l^ one has that 

ei^g{x) = 7r 2 / -e 1^1 dy = tt 2e*l''l / e dy 

Jr3 {x - y) {x - y) 

= 7r-f^/ (l-iA,)e-2-^^ -dy 

{x}^ 4 {x-y) 

= ^"^7Tt/ e-^^'-Ml - X) ^ -]dy 

{x}^ 4: ^ I {x-y) i 

The integral in the last term is bounded uniformly with respect to x. Indeed, all derivatives 
of e^*""*^^!^! are integrable and all derivatives of -^jz^ are uniformly bounded in x and y. 

We deduce that \e^'^g{x)\ < C(x)~^ which implies that e^'^g G L^. This completes the 
proof of the lemma. □ 

We state now our stability result for small solutions. 

Theorem 4.3. There exists an absolute constant £3 > with the following properties: 

• If f,uo G are such that 

(4.5) ||A-ip/||i3,.o + ||no||L3,oo <e3 

then there is a unique solution u G L°^{M.^, L'^'°^) of J^.ip , weakly continuous with 
respect to t £ [0,oo), satisfying 

(4.6) sup \\u{s)\\l3,oo < 2||tio||L3,oo + 4||A"^P/||i3,oo. 

s>0 

• Let p G (|,oo) and suppose in addition to (14. 5 p that uq G If u is the above 
solution then, 

u G L°°(M+, LP'°°) if and only if A^^P/ G L^'^. 
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• Let p G (|,oo), and q > min{3,p}. Suppose in addition to (j4.5p that uq G LP'°° 
and A~^Ff G L^'^^. Let also U G L'^'°° n L^'^^ be the unique stationary solution 
given by Theorem \2.2\ (see also Remark \ 2. 8\) . (We assume here that £3 < ei, the 
constant introduced in Theorem \2.S\) . 

(i) There is a function e{q) > such that if £3 < e{q) then, for some constant 
C>0, 

(4.7) \\u{t) - U\\li < Ct"5(=EW)-i\ Vmin{p, 3} < g < 00. 

In particular, u{t) — [/ — )■ in L'^ as t — )■ 00 for all q > min{3,p}. 

(ii) // I < p < 3, then u{t) U weakly in LP'°° as t — )• 00. Moreover, u{t) — )■ U 
strongly in LP'°° if and only if e^^{uQ — U) ^ in L^'°°. 

(iii) // 1 < p < 3, then the conclusion of the previous item can be strengthened as 
follows: 

(4.8) \\u{t) -U - e*^(uo - U)\\li < Ct^^^^l 

for all < q < p and for some constant C > independent of t. 
In particular, u{t) —[/—)■ in L"? if and only if e^^{uQ — U) ^ in L'^ as t ^ 00, 
for all^<q<p. 

Notice that in (|4.7p neither u{t) nor U belong in general to L'^. Similarly, the terms 
appearing in the LHS of (j4.8p in general do not belong, separately, to L''. In other words, 
the difference u{t) — U is better behaved than the solutions themselves. 

Remark 4.4. In the particular case p = q = 2, the preceding theorem contains an in- 
teresting variant of the stability result for finite-energy solutions obtained in [3] with a 
different method. Indeed, consider a stationary solution [/ G n L^'°^ and a perturba- 
tion wq £ L"^ n L^'°°. According to conclusion (iii), the solution u of the non-stationary 
Navier-Stokes equations starting from uq = U — wq satisfies, under the above smallness 
assumptions, u{t) — )• [/ in as t — )• 00 (we use here that c^^wq — )■ in L^). Explicit 
convergence rates can be given, e.g., if the perturbation belongs to additional function 
spaces. For instance, when wq £ L2'°^ r\ L^°°, then 

\\u(t) - Uy < Ct-^l^, as t ^ 00. 



Remark 4.5. Let us present some further immediate consequences of this theorem. If the 
perturbation satisfies G L^, then e wq — )• in and so in as t — )■ 00. This in 
turn implies, by (ii), 

u{t) -^U in as t 00. 

More generally, according to the second part of Lemma [4.2( such conclusion remains valid 
when lim i?mes{|wol > R}^ = 0. However, notice that neither wr, G L'^'°° is sufficient 

nor lim i?mes{|wo| > R}^ = is necessary to ensure this result. 

3 1 
In the same way, in the case | < p < 3, the condition hm iimes{|wo| > R}^ = 

implies that u{t) U G LP'°° as t — t- 00. But in this case the stronger condition wq G V 
would imply also, by (iii), the stronger conclusion u{t) — )■ [/ in L^. 
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Remark 4.6. The proof of Theorem 14.31 below will show that Equation (j4.8p holds true for 
the wider range max(l, |) < g < -0^. We did not state the full range for q because the 
most interesting case is q < p and also because it would require showing that in the case 
p < 3, the statement (i) is true with a constant e{q) independent of q. This additional fact 
is easy to prove with a recursive argument, but since it is not really necessary we prefer 
to skip it. 

Proof of Theorem \4-3\ We estimate the forcing term in ()4.ip by integrating the heat 
kernel in time then relying on a fixed point argument making use of Lemma 14.11 The 
relation 

(4.9) r e(*-^)^ ds = e*^A"i - 

Jo 

that follows since both operators have the same symbols, gives 

(4.10) II re(*-^)^P/ds||^3,o. < 2||A-ip/||^3..o. 

Jo 

We used above that e*^ is a convolution operator with a function of norm equal to 1. 
Given ()4.10p . the first part of this theorem follows from the work of Cannone and Karch [9]. 
But the proof takes only a few lines, so we give it for the sake of the completeness. 

Using again that the kernel of e*'^ is of norm equal to 1 we deduce that ||e*^uo||L3,(x) < 
||mo||l3.°o- Therefore, if we denote uq = e*^uo + /q e^*~'^^^P/ one has that 

II'"o||l°°(R+;L3.°°) < lko||L3.°° + 2||A"-^P/||2,3,oo. 

To apply the fixed point argument we introduce the notation B{u,v) = Lq{u® v) and 
rewrite (|4.ip as 

(4.11) u = uo-B{u,u) 

The bound (j4.2p . with p = 3, and hence 9 = | combined with the Holder inequality from 
Proposition 11.11 yields 

\\B{u,v){t)\\L3,'^ < C(supC||n(s)||i3,^ ) (sup||t;(s)||2,3.oo ) 

^s>0 ' ^s>0 ' 

We apply this estimate combined with the fixed point argument given in Lemma 12.61 
to the operator $(ti) = mq — B{u,u) in the space L°°(M+,L3'°°). This approach yields 
the existence of a unique solution u € -L°°(M+, L^'°°) provided 4C||no||Loo(]g^.j;^3,oo') < 1. 
Lemma [2.61 also insures that the solution satisfies ||w||j;,oo(]k_^.j;,3,oo) < 2||no||j;,oo(]K_^.j;,3,oo), 
establishing part 1 of the Theorem. 

To prove the second part of the theorem we establish first the cases p £ [2,7] with a 
fixed point argument then treat the other cases with an interpolation argument. First, 
combine ()4.2p with the Holder inequality to establish 

(4.12) \\B{u,u){t)\\LP.oo < C(p)fsup||n(s)||LP.oo)(sup||u(s)||i3,oo) 

^s>0 ^ ^s>0 ^ 

Let C be the maximum value of the constant in the above equation for p € [2, 7], we require 
8^3 < 1/C. Considering again the sequence of approximate solutions (ui) constructed in 
the usual way, and making use of (14. 6p we see 

||'Ui+l(t)||LP,oo < sup ||uo||lp.°° + ^SUp||ni(s)||LP,oo 
s>0 ^ s>0 



LP-SOLUTIONS OF THE STEADY-STATE NAVIER-STOKES EQUATIONS (June 1, 2010) 



21 



From this estimate the "if" statement in the second claim follows for p € [2, 7]. 

If j5 G (|i2), through interpolation we find that for all r G (2,3) we have that uq S 
L°°(R_l_; L^'°°) and therefore u S -L°°(M+; L^'°°). Appealing to (j4.2p and again combining 
it with the Holder inequality we see 

(4.13) \\B{u,u){t)\\LP,o^ < Csup||'u(s)|||.,oo 

where r = G (2,3), hence the right hand side is bounded. Combining this estimate 
with (|4.1ip is enough to prove the "if" statement in the case p G (f ,2). If p G (7, oo) we 
again interpolate to get u G L°°(M+; L'''°°) for all r G (3,6). Choosing again r = G 
(3, 6) in (I4.13P finishes the "if" statement in the second claim. To establish the "only 
if" part of the claim combine (|4.12p with (|4.1ip and notice that the RHS of (|4.9p tends 
to — A"-*^ as t — )■ oo. The weak continuity u{t) — >■ u{t') for t ^ t' and t' G [0, oo) (the 
continuity is actually in the strong topology of L^'°° for t' G (0, oo) ) is proved as in 



It remains to prove the third part of the theorem, the stability results for stationary 
solutions. We begin with Claim (i). Let q > min{3,p}. It is worth noticing that for 
g > 3, a stability result in the L'''°°-norm, as well a decay estimate of the form \\u{t) — 

^^(i)||L9 < Ct 2^ 1 was stated in [91 Proposition 4.3]. However, it seems that the 
argument briefly sketched in cannot be directly applied to the case where the second 
solution v{t) is stationary, because a non-obvious generalization of Lemma l4.ll would be 
needed. Therefore, we provide a detailed proof of estimate (j4.7p . 

Let w = U — u and = U — uq. Then this difference w satisfles the mild PDE 

(4.14) w{t) = e^^w^-f e'^^^'^^^V ■ {u®w + w®U){s)ds. 



Moreover, our smallness assumptions on no and / and the usual fixed point Lemma 
imply that w can be obtained as the limit in L°°{W^ , L'^'°°) of the approximating sequence 
(ujfc), defined by 

Wk+i = e^^uP - B{u, Wk) - B{wk, U), 

where the recursive relation starts with wo{x,t) = e^^w^. Moreover, this sequence {wk) is 
bounded in L°°(M+, LP'°°). 

By the semigroup property (recall that F{x,t) denotes the kernel of e*^Pdiv): 

B{u, v){t) = e^^''^B{u, v){t/2) + f F{t-s)*{u(g> v){s) ds. 

Jt/2 

We deduce, 

Wk+i{t) = e'^w^ - e'^/^B{u, Wk){t/2) - e'^/^B{wk, U){t/2) 

t r-t 

F{t - s)*{u® Wk){s) ds- Fit - s)* (wk U){s) ds. 



It/2 Jt/2 

Now let r = min(3,p) and denote 

M = maxlWu)^ \\ Lr,oa , ||C/||/^r,oo , sup ||u(s)||Lr,oo}. 

s>0 

By Lemma l4. II and using that the sequence Wk is bounded in L°°(M+, L^'"^) 
\\Biu,Wk)it/2)\\Lr,o. + \\Biwk,U)it/2)\\Lr,.. < CrM. 
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A heat kernel estimate now implies, for all g > r and for some constant > independent 
of q, 



\Wk+l 



LI'- 



+ 



t/2 



F{t — s) * (n 05 Wk){s) ds 



Li'' 



+ 



t/2 



F{t - s)* {wk (S) U){s) ds 



Li'' 



From Lemma |4. II with Holder's inequality we have 



t/2 



F{t — s) * {u^ Wk){s) ds 



+ 



t/2 



F{t - s) * {wk(S)U){s)ds 



LI'- 



<Cqe3 sup \\wk{s)\\Li 

s&[t/2,t] 



Let 



then 



Wk{t) 



sup \\Wk{T)\\Li,< 
T£[t,oo) 



Wk+i{t) < C'.Mt 
Iterating this inequality implies 

fc-i 



l^r q' 



n=0 



3/1 3/1 



-{---) 



2^r q' ] 



t 2 



3(1„1) 

\ J- q J 



provided 



C^3 2^ 



q' < 



3(l_ll 1 

A slightly more stringent smallness condition and independent on r > | is, e.g.. 



(4.15) 



£3 < e(g) := 



Now assuming ()4.15p and letting A; — t- oo we get, 



(4.16) 



\w 



LI'' 



< 2C.j.Mt ^^-r i> ^ for r = min(3,p) and q> r. 



Writing the above estimate for q — rj and q + i], for some r] > small enough and interpo- 
lating the L'^-space between L'^~'^'°° and L'^'^'^'°° shows that the above estimate remains 
valid with ||ii;(t)||2,9 on the left-hand side. This establishes the stability result (j4.7p . 

We now prove Claim (ii). The weak convergence u{t) ^ [/ in for | < p < 

3 is obvious since the solution u[t) is bounded in LP'°° and goes to U in the sense of 
distributions (even in L*^, g > 3, as implied by the previous part of the proof). 

On the other hand, the proof of the necessary and sufficient condition for the strong 
convergence result in the L'^'°°-norm is given in [9l Theorem 2.2] and in [9l Corollary 4.1], 
hence we will skip it. The necessary and sufficient condition for the strong convergence 
result in the LP'°°-norm, with | < p < 3 is a direct consequence of Claim (iii) which we 
now prove. 
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Let now | < p < 3 and ^ q ^ P- Given these restrictions, there exists some 

qi G [p, 4] such that the fohowing relations hold true: 

1111 1 ,1111, 
< — < -, — < min(l - -,- + ). 

q p qi p qi ^ p 3 q p' 

We go back to the equation for w given in (|4.14p . We estimate — wo{t)\\Lq using 

Propositions 11.11 and \1.2\ the bound ||tL'(t)||L9i.°° < Ct n' (consequence of (I4.7p with 

q = qi) and the the fact that U G 17'°° and u £ L°°((0, oo), LP'°°). We get 

p n'\\{u0w + w ®U){s)\\ _qip_^ds 
1 



w{t)-WomL'> <c 


/'(.- 






< c 


/o 






< c 


'0 




l3(1 



'p n^||w(s)||L«i,-(||n(s)||LP,- + \\U\\ ) ds 

.1_J_) _3cl_J_\ 
P 91 '5 2Vp 11' ds 

1 I 3 _3 
< (7^2 + 29 p. 

The theorem is now completely proved. 

We finally show our stability result for large solutions. In lay words the next theorem 
shows that global solutions u of the non-stationary Navier-Stokes equations with a small 
constant forcing term / tend to forget their initial data. Indeed, such solutions will 
converge to the small stationary solution U driven by the same forcing / even if they are 
initially large. 

Theorem 4.7. There exists an absolute constant £4 > with the following property. Let 
u G -L^^([0, 00); L^'°°) n -L^^^([0, 00); L^) he a global solution of the evolutionary Navier- 
Stokes equations with a constant in time forcing f such that A~^P/ G L^'°° H L'^ and 

(4.17) A(uQ,f) = limsupiimes{|uo| > i?}^ + || A"^P/||i3,oo < £4. 

R-^O 

Let U G L^'°°nL'* he the unique stationary solution constructed in Theorem \2.2l Then we 
have that 



limsup ||n(t)||j;^3,oo < 22A{uq, f); 



00; 



u[t) U weakly in L^'°° as t 
u{t) U in L'^'°° as t —7- 00 if and only if e*^(no — C/) — t- strongly in L^'°° 
t — )• 00. 



as 



Proof. The idea of the proof is the same as in [15] where it was proved that any global 
solution of the Navier-Stokes equations without external force goes to in the Besov 

spaces Bp^q ^ when the time becomes large (see also [l] for the case of VMO'^^). It 
consists in decomposing the initial velocity in a small part plus a square integrable part. 
The small part remains small by the small data theory and the square-integrable part will 
become small at some point by using some energy estimates. 

Here we use Lemma l42] to decompose uq = vq+wq where vq G L^nL^'°° and ||t(;o||L3,oo < 
2A{uo,f). Assuming that 3£4 < £3 where £3 is the constant from Theorem 14.31 we can 
apply that theorem to construct a global solution w of the Navier-Stokes equations with 
forcing term /, initial velocity wq and such that 

|h(t)||L3.oo < 8A{uo,f) for ah t > 0. 
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Moreover, according to relation (j4.7p the solution w satisfies the following decay estimate 
supts — U\\n < oo. Since f/ G we infer that w G Lf^^{[0,oo);L'^). 

The difference v = u — w verifies the following PDE: 

(4.18) dtv- Av + u-S/v + v-Vw + Vp =0 
whose integral form reads 

(4.19) v{t) = e^^vo-[ e^^-'^^^FV ■ {u(g>v + v(g>w){s)ds. 



We show first that v G C''^([0, oo); L^). The first term on the RHS above clearly belongs 
to this space. We show that so does the second term. The kernel F{t) of the operator 
e*^Pdiv is of the form F{x,t) = ^^^(^,1) with F(-,l) G n C LP'^ for aU 

1 < p < oo and 1 < q < oo. In particular, ||-F(t)|| 6 2 < t~i so that F G Lj^^{[0,oo); L^'"^). 

3 

By the Holder inequality we also have that u^v + v0w(^ L^^([0, 00); L2 '°°). Since the 
last term in (j4.19p is the space-time convolution of F with u v + v w, we infer that it 
belongs to C°([0,oo);L2). 

For < 5 < 1, let J5 be a smoothing operator that multiplies in the frequency space 
by a cut-off function bounded by 1 which is a smoothed out version of the characteristic 
function of the annulus {6 < < ^}. We also introduce an approximation of the identity 
ip^ in time. 

Given the additional regularity found for v above, we remark that we can multiply the 
equation of v expressed in (|4.18p by * * JgV and integrate in space and time from to 
to t, with tQ > 0, to obtain that 

(4.20) \\ip^ * J5v{t)\\l2 + 2 I * Jsv{s)\\l2 ds = * Jsv{to)\\l2 

Jto 

+ 2 / U ■ ^{fri * (frj * Jsv) ■ V + 2 / V ■ V {ifj^ * 'Prj * ^5^) ■ W. 

Jo Jr^ Jo jr3 

We let now r/ — )• 0. Given the time continuity of v with values in L^, we have that 
iprf * J5v{t) — )• Jsv{t) and tprj * Jsv(to) — Jsvito) in as 77 — )• 0. The other terms in (|4.2Up 
pass easily to the limit r/ — t- 0. Therefore, taking first the limit 77 —t- in (|4.20p . and second 
to ^ and using again that v G C^[0, 00); L^) we get that 

(4.21) \\Jsv{t)\\l2+2 [ \\VJsv{s)\\l2ds=\\Jsvo\\l2 + 2 [ [ u-VJjvv 

Jo Jo 

+ 2 / vVJjvw 
Jo Jr^ 

We bound the last two terms on the RHS as follows 

2 /" /" u-Vjjv-v + 2 [ [ vVJsV-w<2[ \\V Jsv\\l^\\v\\l4\\u\\l^ + WwWl"^) 
Jo Jo Jo 

<\ I \\^Jsv\\l2+ [ Ibiii4(i|w|ii4 + ii^/;|ii4). 

^ Jo Jo 

Plugging this in (|4.2ip yields 

\\J5v{t)\\h + r \NJ5v{s)\\l2 ds < WJsvoWh + r ||^||i4(||n||i4 + \\w\\l,). 
Jo Jo 
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Since u,v,w G Lf^^{[0,oo); L"^), the RHS above is uniformly bounded with respect to 6. 
Letting 6^0 imphes thanks to the Beppo-Levi theorem that ||Vf(s)||^2 ds < oo, that 
■isv€Ll^{[0,ooy,H^). 

We go back to (14.2ip and estimate 

(4.22) 2 [ [ u- Vjjv ■v = 2 [ [ u- V j|w • (1 - JI)v 
Jo Jr^ Jo Jr^ 

< 2||u||i4(o,t;L4)||V7;||i2(o,t;L2)||(l - Js)v\\L'i{o,t;L-i) ^0 as (5 ^ 0. 

We observe now that if^(]R^) ^ L^'^(]R^). This imbedding follows from the Young 
inequality for Lorentz spaces after noticing that (— A)~2 is a convolution operator with a 
function bounded by which therefore belongs to L2'°°. We use this fact together with 
the Holder inequality to bound the last term in (j4.2ip as follows 



2 

L2 



(4.23) 2 /" /" vVjjv-w<cl ||-i;||i6.2||Vj|u||i2||it;||i3,oo < C^(no,/) /" ||Vt;| 
Jo Jr? Jo Jo 

Using (|4.22p and (|4.23p in (j4.2ip . letting 5 ^ Q and using the Beppo-Levi theorem we 
infer that 

rt 



\v{t)\\l, + 2[ \\Vv{s)\\l,ds<\\vo\\l,+CA{uoJ) [ \\Vr 
Jo Jo 



|2 

Il2 



If we further assume that Ce4 < 1, then CA(uq, /) < 1 so the relation above implies that 
V £ L°°(M+;L^) n L'^{R+; H'^). By interpolation and from the imbedding C L^'°° 
we infer that v £ L^(R+;H^) C L^(M+; L^'"^). So there exists a time T = T(e4) such 
that ||f (r)||^3,cx) < A{uo,f). Since we also have that ||u)(r)||^3,oo < 8^4(^0,/) we infer 
that ||n(T)||j;^3,c3o < 9A{uo,f). Assuming further that 10^4 < £3, Theorem 14.31 allows 
to construct a small solution starting from time T, a solution whose norm will be 

bounded by 22^(uo,/). We will prove below a uniqueness result stating that u must be 
equal to this small solution starting from time T. Once this is proved, the first part of 
the theorem follows. Moreover, using again that our solution u becomes small after the 
time r, the second and the third part of the theorem are consequences of Theorem 14. 3[ 
Except that the equivalent condition for u{t) to converge strongly to U in L^'°° is that 
e^^{u{T) — ?7) — )• strongly in as t — )■ 00. To finish the proof it therefore suffices to 
show that 

e*^(n(r)-C/)*-^0inL3'°° ^ e*^(no - C/) *^ in L^'"". 
This is a consequence of the following sequence of equivalence relations: 

e*^(n(r)-?7)*^0inL3'°° ^ e*^(ii;(r) - [/) in L^'^ 

^ e*^(txo-C/)*^OinL3'-. 

We used above that v{T), vq G and the decay estimates for the heat equation to deduce 
the first and fourth lines of the relation above, and Theorem 14. 31 twice for w, starting from 
time t = and from time t = T to deduce the second and third lines. This completes the 
proof provided that we prove the announced uniqueness result. 
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Let u be the small solution starting from time T with initial velocity u{T) constructed in 
Theorem l4.3l and set v = u—w. As above we have that v G C'^([T; oo); L^)nL^([T; oo); H^). 
Then v — v solves the following equation: 

dt{v — U) — A(u — v)+u- V{v — v) + {v — v)- Vu = — Vpi. 

As in the previous argument, one can prove that this relation can be multiplied hy v — v 
and integrated from T to t to get that, for all t >T, 

\\{v - v){t)\\l, + 2 £ \\V{v - v)\\l, = £ l(v-v).V{v-v)-u 

<C f\\V{v-v)\\l4u\\L^,o.<CA{uo,f) f\\V{v-v)\\l,< f\\V{v-v)\\l, 

JT JT JT 

provided that A{uo,f) is sufficiently small. We infer that v{t) = v{t), that is u{t) = u{t) 
for all t > T. This completes the proof of the theorem. □ 



Remark 4.8. We also have stability in U''°° for large solutions. More precisely, suppose 
that in addition to the hypothesis of Theorem 14 . 71 we assume thattto G LP'°° withp G (|,3). 
Then u £ L°°(M+; LP'~), u{t) U weakly in and u{t) - e*^(«o - U) ^ U in LP^°° 
as t — 7- oo. This follows easily after applying Theorem 14.31 starting from the time T when 
the solution becomes small. One only needs to show the following two facts: 

• if ito G and A'^P/ G LP'°° then u G L°°(0,T; LP'°°); 

• e(*-'^)^(?x(r) - U)- u^^{uo - [/) strongly in LP'°° as i ^ oo. 

To prove the first assertion, we observe that, with the notation from the proof of Theo- 
rem [43] (namely the notation used in relation (14. lip ) one has that uq G L°°(]R_|.; L^'°°). 
Moreover, by the Holder inequality and using the standard decay estimates for the heat 
equation we can bound 



(4.24) \\B{u,u){t)\\LP < [ We'^^-'^^fV ■ {u^u){s)\\Lpds 

Jo 



<C I (t-s) '^^\\u{s)\\l3.^ <Ct sup ||ii(s) 11^3.00. 



10 0<s<t 

We infer that B{u,u) G L°°(0,T;LP) C L°°(0, T; LP'°°), so by we also have that 

u G L°°(0, T; LP'°°). To show the second assertion, we observe that it is sufficient to prove 
that 

n(T)-C/-e*^(uo-C/) GL9'°° 
for some q < p. But u — U verifies the PDE 

dt{u - U) - A{u - U) + u ■ Vu - U ■ VU = -Vp2 
whose mild formulation implies that 

u{T) -U-e^^{uo-U) = - [ e(^~")^PV ■ {u u - U U){s) ds 

Jo 

The same estimate as in (j4.24p shows now that the RHS belongs to L'^'°° for any | < g' < 3, 
in particular for some q < p- 

Moreover, if no G L^, then the previous argument shows that u G L°°{0,T; L^). Prom 
Theorem 14.31 applied starting from time T we infer that u G L°°(M+;LP) and u(t) U in 
LP as t ^ oo. 
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Remark 4.9. We observe that the condition imposed on the initial velocity by the hypoth- 
esis of Theorem 14.71 does not imply that uq is close in L^'°^ to the smooth functions in 
5(M^). Indeed, that would require to have that the quantity limsupi2mes{|uo| > R}^ is 

R^oo 

small too. This condition is not necessary in Theorem 14.71 
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